Tensor interaction constraints from (3 decay recoil spin asymmetry of trapped atoms 
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We have measured the angular distribution of recoihng daughter nuclei emitted from the Gamow- 
Teller (3 decay of spin-polarized *°Rb. The asymmetry of this distribution vanishes to lowest order 
in the Standard Model (SM) in pure Gamow- Teller decays, producing an observable very sensitive 
to new interactions. We measure the non-SM contribution to the asymmetry to be At= 0.015 
± 0.029 (stat) ± 0.019 (syst), consistent with the SM prediction. We constrain higher-order SM 
corrections using the measured momentum dependence of the asymmetry, and their remaining un- 
certainty dominates the systematic error. Future progress in determining the weak magnetism term 
theoretically or experimentally would reduce the final errors. We describe the resulting constraints 
on fundamental 4-Fermi tensor interactions. 

PACS numbers: 23.40.-s,32.80.Pj,14.80.-j 
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I. INTRODUCTION 
A. Search for tensor interactions 

Effective 4-Fermi contact interactions contributing to 
beta decay can be classified by the Lorentz transforma- 
tion properties of the contributing lepton and hadron cur- 
rents [H. The Standard Model (SM) of particle physics 
contains vector (V) and axial vector (A) interactions with 
sign 'V-A'. Extensions to the standard model can pro- 
duce effective scalar and tensor interactions. The ob- 
servable measured here was developed shortly after the 
discovery of parity violation by Treiman, who realized its 
sensitivity to tensor interactions [2]. 

To lowest order and neglecting the Fermi function, the 
angular distribution W[0\ of the daughter nuclear recoils 
with respect to the nuclear spin, integrated over all final 
momenta, is given by , 
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with / the total nuclear spin and M the spin projection 
along the quantization axis, so that P is the vector po- 
larization of the nucleus, and T is the 2nd-rank tensor 
alignment. So the recoiling daughter nuclei detected in 
singles from the P decay of polarized nuclei have spin 
asymmetry, integrated over all recoil momentum, Aj-ocoii 
= -{Ai3+B,j) where Afj and B^, are the f3 and 1/ asymme- 
tries. The alignment coefficient c depends on the nuclear 



spins and is calculated in Ref. Q and Section ITC] 

The coefficients xi and X2 are calculated from integra- 
tions over the momentum of the other outgoing particles, 
and their detailed dependence on the energy release Q is 
given in Appendix A. In the limit of high Q, as is the 

case for our *°Rb decay, xi 5/3 and X2 1/2. 

(Note that xi vanishes as Q ^ 0, as it must because the 
helicity of the /?"*" vanishes as its momentum decreases to 
zero.) 

For Gamow- Teller decays in the SM, Ai3=—B^, so 
^rccoii=0 to lowest order. Vector and axial vector cur- 
rents that couple to right-handed neutrinos also cancel in 
the sum Ap+Bi, [1]. That makes Arecoii sensitive only to 
fundamental lepton-quark tensor interactions: it is nei- 
ther sensitive to effective scalars nor to right-handed cur- 
rents. The exclusive sensitivity to tensor interactions is 
only true for one other observable, the 13-v correlation in 
pure Gamow- Teller decay [1, @ . 

The near-zero value of Arccoii in the allowed approxi- 
mation for pure Gamow- Teller decays also makes it a very 
attractive experimental observable; for example, the po- 
larization P does not have to be known to high precision 
to extract the new physics. We will see below that recoil- 
order corrections produce effects ^ 0.01 in the absence of 
new physics. Note from Appendix A that the recoil spin 
asymmetry remains proportional to Afj+B^, even before 
the integration over recoil momenta. 



1. Explicit sensitivity to 4-Eermi tensor coupling 
coefficients 

The tensor-axial vector Fierz interference term is lin- 
ear in a particular combination of the 4-Fermi tensor cou- 
pling coefficients Ct and [3j: 
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This term appears in the full expression for the decay dis- 
tribution before integration over the /3 momentum [3| . It 
is well-constrained by other experiments in lower-energy 
/3 decays [7], so we will cite results below both with it 
unconstrained and with it constrained to be small. Its 
inclusion complicates the expressions, so we include it 
in our full simulations by numerical integration; we also 
show the full expression (without Fermi function) in Ap- 
pendix A. Note that Eq. 1 has consciously assumed the 
Fierz term to be zero [4I to simplify the expressions. 

In absence of this Fierz term, the contribution to Arccoii 
from non-SM interactions becomes just the product of 
the 4- Fermi tensor constants Ct and j^]: 
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where the ± sign is for j3'^ decay. The coefficient A/j' is 
given in Ref. [3] (see Section ITC|l . 

The combination Ct + C'^ describes a tensor interac- 
tion that couples to standard model left-handed neutri- 
nos, while the coupling Ct — C'j, describes an interaction 
coupling to neutrinos with non-SM helicity Q. So the 
sensitivity of Arocoii to CtO^ produces sensitivity to in- 
teractions with both SM and non-SM chirality. 



B. Tensor interactions in other experiments and 
theory 

Until recently, individual nuclear /5-decay correlation 
experiments and global fits have been consistent with the 
standard model without a tensor interaction [13, [lH . 
A recent global fit of nuclear and neutron (3 decay 
data including scalar and tensor terms coupling only to 
standard-model left-handed v's (i.e. assuming Ct=C^) 
gives Ct/Ca = 0.0086±0.0031, while excluding the lat- 
est neutron lifetime measurement [l^ from that fit brings 
the result into agreement with the Standard Model at 
one standard deviation Q. Such difficulties in combin- 
ing many experiments with different systematic errors 
in global fits can be avoided by dedicated experiments 
sensitive only to tensor interactions, like the present ex- 
periment. So the present measurement becomes useful if 
it achieves '--^ 0.01 accuracy in the recoil asymmetry. 

A tensor interaction coupling to right-handed neutri- 
nos would produce a contribution to the mass of the 
standard model left-handed neutrinos. An order-of- 
magnitude calculation suggests that if \Ct — C!p\ were 
~ 0.02, such an interaction would account for neutrino 
masses ^ 3 eV [l3| . the present experimental upper limit. 
This provides a motivation for direct correlation measure- 
ments in order to constrain this possible contribution to 
neutrino masses. 



The PIBETA collaboration reported a statistically sig- 
nificant deviation from the standard model [14] in Tr—nyej 
decay that could be explained by a finite tensor interac- 
tion. The same group has made further dedicated ex- 
periments and now eliminated the possibility of a tensor 
interaction 15]. The most restrictive limits from ir^vej 
decay are on interactions with same chirality as the Stan- 
dard Model, i.e. constraining a Fierz interference term 
that is linear in the small tensor term, although the op- 
posite chirality is also considered in Refs. [1, [la]. Sensi- 
tivity in the present experiment ^ 0.001 level would be 
needed to reach the sensitivity probed by the tt decay 
experiments. 

A renormalizable tensor interaction can be generated 
by the exchange of spin-0 leptoquarks [13] • In an ex- 
ample of an explicit model, a recent analysis of the possi- 
ble one-loop corrections in SUSY models shows they can 
produce tensor and scalar interactions as large as 0.001 
in the Fierz interference term 6^ from left-right sfermion 
mixing in the first generation, physics that is otherwise 
difficult to constrain "ISl. 



C. Rb and recoil-order corrections 

The decay of F^=l+ *"Rb is primarily to two states, 
74% to the ground 0+ state and 22% to the first excited 
2+ state (see Figure [T]). The coefficient Xij' in Eq. [5] 
is 1 for the l+^O"*" transitions and -1/2 for the 1+^2+ 
transitions. So when averaged over the transitions, the 
average Xjj/ for ^°Rb is 0.64. The coefficient c in Eq. 1 
is 1 for the l+^^O"*" transitions, and -1-0.1 for 1+— >2+. 

The 1+— >0^ transition has higher-order corrections in 
the SM proportional to the weak magnetism form fac- 
tor 6m and the induced tensor form factor d [l^. The 
contributions to the angular correlation coefficients from 
these recoil-order terms are scaled by Efj/Maudcon, and 
they produce a small nonzero recoil asymmetry within 
the Standard Model. 

These recoil-order terms are given by sums over indi- 
vidual nucleon matrix elements 



bM/A = smMgt + {f\\J2 ^k^k IN) (6) 

k 

d/A ^ gA{f\\J2r+zdkXlk\\i) (7) 

k 

using the notation of Ref. ^2d\ , which has detailed consid- 
eration of the recoil-order matrix elements needed here. 

The first term in the expression for weak magnetism, 
Lm/A, is from the anomalous isovector magnetic moment 
of the nucleons, and is given by the dimensionless num- 
ber gM—4:.7 multiplying the Gamow- Teller matrix ele- 
ment Mqt- The second term in Bm/A requires a detailed 
nuclear structure calculation with a large fractional un- 
certainty in this mass region, as deformation effects make 
shell model calculations difficult. Weak magnetism terms 
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II. EXPERIMENTAL TECHNIQUES 

We describe below the trap apparatus which we use to 
polarize the ^''Rb nuclei and measure the angular distri- 
bution and momentum of the ^''Kr daughter. 

The collection and trapping of the ^"Rb used a two- 
very similar to that in previous /3- 
. The ^''Rb 30 keV ion beam from 
the IS AC facility at TRIUMF was w 2 x lOVsec from 
a zirconium carbide target. The ion beam was stopped 
in a neutralizing foil made of zirconium. The resulting 
atoms were collected in the first magneto-optical trap 
(MOT) [l^l . The trapped atoms were transferred to the 
second MOT to minimize backgrounds [2l| and provide 
an environment for polarization. The number of atoms 
continuously trapped in the detection trap was « 2 x 
10^. 



FIG. 1: The /3 -decay scheme for Rb, showing literature 
values of the spin, parity, branching ratio, logio(ft), and ex- 
citation energy of the final levels in ®°Kr, along with the Q- 
value (maximum kinetic energy) for (3'^ decay and the parent 
half-Ufe. Energies in keV. 



for a variety of Gamow- Teller decays have been calcu- 
lated and measured, and the values generally do not de- 
part greatly from the value of 6a^/(^Mgt)=4.7 [HllOl. 

The induced tensor term d is more poorly character- 
ized. We will see below that it contributes values to the 
experimental recoil asymmetry that are roughly constant 
with recoil momentum, so we will be able to fit for d si- 
multaneously with the non-SM tensor interaction. 



D. Transitions to excited states 

The 22% 1+ to 2+ transition in principle adds some 
complication to the recoil-order matrix elements. The 
first-order form factors, bm and d, can be different from 
those for the transition to the ground state. In addition, a 
total of five additional form factors appear at 2nd-order 
in recoil terms /Mnucicon • In this version of the ex- 
periment, we will let bM and d float phenomenologically 
below, so that we implicitly include b^i, d for the excited 
state in our analysis. In other words, bM and d for both 
the ground and excited states are treated as producing 
the same functional dependence of the recoil asymmetry 
on recoil momentum. This is an excellent approximation 
as the effects on the recoil asymmetry have very similar 
momentum dependence compared to the ground state, 
particularly once they are averaged over the momentum 
spread induced by the final 7-ray emission. The values 
of 5m and d that are extracted are then weighted aver- 
ages of the transitions. We ignore the 2nd-order terms in 
E^/Mnucieon, which is a good approximation at the level 
of accuracy reached in this version of the experiment. 



A. Detection geometry 

The detection geometry is shown in Fig. [51 The time- 
of-flight (TOF) of the daughter nuclei from nuclear f3 
decay in singles (i.e. not in coincidence with the /?) can 
be measured by using the atomic shakeoff electrons as 
a trigger, as developed by LBL researchers [1^. A uni- 
form electric field of average value 800 V/cm collects ions 
produced in f3 decay to a 25 mm diameter microchan- 
nel plate (MCP) Z-stack for time readout, backed by a 
position-sensitive resistive anode. The electric field also 
collects the atomic electrons to a second MCP detector 
on the opposite side. A positive ion produced in /3+ de- 
cay will 'shake off' at least two atomic electrons, and 
electrons up to ~100 eV energy are completely collected 
by the field into the detector. So this technique increases 
the efficiency for recoil detection by a factor of about 30 
compared to /J"*" detection (which had solid angle 1% in 
our geometry [2?|). This produced approximately 100 Hz 
of daughter recoils in coincidence with the electron de- 
tector, for the average number of atoms trapped of about 
2x10^. 

The electron detection geometry and electric fields are 
optimized to ensure efficient detection independent of the 
electron's initial kinetic energy and angle. We measured 
« 35% detection efficiency for low-energy electrons from 
the laser photoionization of ®°Rb (Fig. [3]). By photoion- 
izmg stable **^Rb, we have reproduced data in the lit- 
erature showing that the electron efficiency of MCPs is 
optimized at 500 eV impact energy, and changes by less 
than 10% between 500 and 1000 eV |28i]. The shakeoff 
electrons are expected to have kinetic energies similar to 
their atomic binding energies, a few lO's of eV f29|, so we 
arrange the fields so that the electrons impact the MCP 
with approximately 500 eV more than their original ki- 
netic energy. To make the electric field more uniform in 
the region traversed by the electrons than in the geome- 
try of Ref. 22], the grid in front of the electron detector 
is biased to a potential close to that of the final electrode 
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FIG. 2: (Color online) TRINAT detection apparatus. The 
MCP for ion detection has position-sensitive resitive anode 
readout. Given the known uniform electric field, measure- 
ment of the recoil ion TOF and impact position determines 
the initial momentum. The apparatus is similar to that 
of [12], with an MCP for electron detection added for the 
^''Rb decay asymmetry. The *°Rb is spin-polarized by op- 
tical pumping with light resonant to the Dl transition from 
a 50 mW diode laser. The polarization is monitored with (3 
AE — E phoswiches using plastic/CaF2 (Eu) fast /slow scin- 
tillators (which are at -30 and 150 degrees out of the plane 
shown here). 



in the field assembly. 
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FIG. 3: (Color online) Energy diagram of the atomic levels 
for optical pumping of the Rb atom, along with the pho- 
toionization scheme for position determination. The atomic 
hyperfine structure is not to scale. The frequency splittings 
were taken from [2711. 



B. Polarization techniques 



The ^"Rb atoms are polarized by switching off the 
MOT light and optically pumping with light at the Dl 
transition (see Fig. [3]) for 30 ^s. Then the MOT light is 
switched on again for 30 /is to keep the atom cloud from 
expanding. The MOT 3 G/cm (horizontal) quadrupole 
field stays on at all times. In a typical MOT with beams 
carefully balanced in power, the atom cloud would be 
centered at zero magnetic field, but such a cloud would 
sample nonzero fields and the polarization would be dis- 
turbed by Larmor precession. In order to add a constant 
magnetic field along the optical pumping axis, we at- 
tentuate two of the beams in the MOT horizontal plane, 
perpendicular to the quadrupole field anti-Helmholtz coil 
axis. The cloud equilibrium position is then at finite B 
field, i.e. at a location where Zeeman shifts produce equal 
absorption from the unbalanced beams. Then we in ad- 
dition apply at all times a uniform field constant of 2.5 G 
with Hclmholtz coils. The result is a cloud equilibrium 
position at the center of the apparatus, with on average 
a 2.5 G field along the optical pumping axis. The cloud 
spatial FWHM of 3 mm samples ±0.5 Gauss of changing 
field. 

Atoms were transferred from the 1st trap every 1.5 



seconds. After each transfer, the polarization state was 
flipped by changing the handedness of the optical pump- 
ing light with a liquid crystal variable retarder. 

The 30 /iS trap on/off polarization off/on duty cycle 
was chosen to minimize the motion of the cloud during 
the optical pumping time, which was a dominant sys- 
tematic in V asymmetry measurements in ^^K [2^. The 
2.5 G bias B field means that there were slightly differ- 
ent Zeeman shifts for the two different polarizations. To 
minimize differences in the atom cloud position, the op- 
tical pumping laser frequency was shifted (by less than a 
linewidth) between the two polarizations. 

The atom cloud position was monitored by photoion- 
izing a small fraction of the atoms with a pulsed laser 
(see Fig. [S]) and recording MCP position and TOF. The 
average trap location was found to shift by 0.030±0.003 
mm with spin flip, which would produce a false asymme- 
try Arecoii of 0.0012±0.0001. The correction was made to 
the data by using the measured trap location when de- 
termining the experimental angle of emission below (Sec- 
tion |TTT|. 
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C. TOF spectrum 

We show a typical TOF spectrum in Fig. 21 decon- 
structed into its components. The different charge states 
are separated in TOF by a uniform 0.8 kV/cm electric 
field. Charge states 1, 2 and most of charge state 3 are 
relatively clean of background. Ions from the 1.4% elec- 
tron capture branch have a large spin asymmetry and 
contaminate the higher charge states, and because their 
asymmetry is large they become a useful probe of the 
polarization (see below). Note that only about 15% of 
the /3+ decays produce positive ions, as opposed to most 
of the electron capture decays. 

We measured the background from /3+'s striking the 
electron detector by lowering the bias voltage of the de- 
tector to exclude the atomic electrons, while keeping the 
grid in front of the detector at the same voltage to keep 
the electric field for ion collection the same. That back- 
ground has the expected large spin asymmetry from (3- 
recoil coincidences, which in this geometry is determined 
by the v asymmetry . The recoil asymmetries shown 
below have been corrected for this background, which 
produces 1.2% of the charge state 1 recoils, 0.59% of the 
charge state 2 recoils, and 0.51% of the charge state 3 
recoils. The average asymmetry correction can be seen 
below in Section [1111 Fig. [HI and will be discussed there. 

There is in addition a 1% background between 1.2 and 
1.4 /is that shows a definite localization on the lower part 
of the ion MCP. This background also appears in natural 
backgrounds and 7-ray source measurements, and may be 
due to an electronic artifact. Although its origin is not 
fully understood, we have measured its spin asymmetry 
to be negligibly small and consistent with zero, so via 
this technique we can correct for its presence to sufficient 
accuracy. 
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FIG. 4: (Color online) (a) Time-of-fiiglit (TOF) spectra for 
recoil coincidences with shakeoff atomic electrons, showing 
decomposition into decay and backgrounds. The small 
electron capture (EC) branch produces large corrections for 
the higher charge states, and is modelled here assuming charge 
state distributions from x-ray photoionization 26]. A back- 
ground of order 1% from /3"'"'s striking the electron detector 
is determined by biasing the detector to exclude low-energy 
electrons (see text). The simulation for charge state -1-3 is 
highlighted; the TOF for -1-3 was cut above 0.63 /is to help 
exclude EC events, (b) Expanded TOF scale of the top figure, 
showing the charge states 9, 10, and 11 that are dominated 
by EC and used for polarization determination (Section lllDfl . 



D. Vector polarization determination 

Here we describe the determination of the nuclear vec- 
tor polarization achieved, 0.55± 0.04. The precision is 
more than adequate because the present observable van- 
ishes, though the final error on the new physics parame- 
ters At and 6t is compromised because it scales inversely 
with the absolute polarization achieved. 

The polarization was optimized by measuring the time 
dependence of the atomic excited state population dur- 
ing the optical pumping, monitored by non-resonant pho- 
toionization with a small pulsed laser (Fig. ^ . The ex- 
cited state population decreases as the polarization in- 
creases: if the atoms become fully polarized, then the 
atom can absorb no more light in this transition and the 
excited state population would vanish. Equilibrium po- 
larization is reached during the last 20 \i& of the optical 
pumping, during which ^recoil is measured. Unlike in 
our previous work with the MOT magnetic quadrupole 
field turned off [1^, we found that the atomic measure- 
ment of the polarization was difficult to quantify with the 



MOT quadrupole field left on. So the polarization was 
measured by nuclear observables. 

The Z?"*" asymmetry was measured using 
plastic/CaF2(Eu) phoswich detectors, in coincidence 
with shakeoff electrons to minimize sensitivity to decays 
from untrapped atoms. The phoswiches are located at 
-30 and 150 degrees with respect to the polarization 
direction, out of the plane of Fig. [51 

The /3 asymmetry for the l"*" to 0+ transition is -fl, 
while for the 1+ to 2+ it is -1/2. So the asymmetry 
grows at the higher /3 energies, as the contribution of 
the 1+ to 2"*" becomes proportionately smaller. Figure [5l 
shows the fit asymmetry as a function of /? momentum. 
For betas in coincidence with recoils in charge states -1-4 
to -1-9, the full solid angle of the recoils is detected, so 
this physical observable is the same as the singles /3 spin 
asymmetry. It has the advantage of being completely 
clean of background from decaying atoms that are not in 
the trap. The result is P=0.53 ± 0.03. 

An entirely different nuclear observable serves as an 
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FIG. 5: Nuclear polarization determined from the (3 asym- 
metry in coincidence with Kr recoils charge states 4-9 in the 
CaF2(Eu) detectors. 



E. Tensor alignment 

We have no direct observables that are very sensitive 
to the tensor-order polarization alignment T (defined in 
Eq. [3]), but we can make adequate indirect constraints 
from the measured vector polarization P=0.55±0.04. For 
nuclear spin I—l and P<1, the population that is not in 
spin projection m/=l must obviously either be in mi— 
or -1. If it were all in mj=-l, then the value of T=-l, 
unchanged from its value for perfect P=l. For imper- 
fect optical pumping spoiled by Larmor precession in the 
quadrupole field, it is much more likely for most of the 
population to be in mi—0; doing this produces our best 
estimate of T=0.35. In the fits, these two extremes do 
not significantly perturb the extracted asymmetry coeffi- 
cients Ai, so we do not mention them further, and simply 
take T=0.35 in the remaining analysis. 




cosfl with respect to polarization 



FIG. 6: Nuclear polarization determined from recoils from 
electron capture (EC) decay (see text), which has a large spin 
asymmetry, (a) Counts as a function of cos 9 for opposite 
signs of polarization, (b) The resulting asymmetry of events 
from (a); the line is a fit with P^i5C=0.29±0.02 (see text). 



additional measurement of the nuclear polarization. Re- 
coils produced in electron capture (EC) decay dominate 
the higher charge states 9, 10, and 11 (see FigureSKb).) A 
further cut on recoil momentum removes a 25% contribu- 
tion from /J^-produced recoils, leaving a clean sample of 
the highest-momentum EC-produced recoils. Figure[6Kb) 
shows the angle dependence (as constructed from the 
MCP position information and time-of-flight) of the re- 
sulting recoils from EC. If the polarization P were unity, 
the EC recoils would have asymmetry unity for the 1^ 
to 0+ transition, and -0.5 for the 1+ to 2+ transition. A 
simple linear fit to the asymmetry in Figure Elb) (there 
is no cos'^e term 0), extracts PAec=0.29±0.02, which 
implies nuclear polarization P=0.57±0.04. 

In summary, the nuclear vector polarization achieved 
from these consistent observables is 0.55 ± 0.04. The 
uncertainty produces an error on the extracted j4i.ccoii of 
7% of its value, which as we will see below is negligibly 
small in this version of the experiment. 



III. ANALYSIS OF RECOIL ASYMMETRY AS 
A FUNCTION OF RECOIL MOMENTUM TO 
EXTRACT RECOIL-ORDER TERMS 

From the ion MCP hit position, ion TOF, trap cloud 
location from photoionization, known uniform electric 
field, and known charge states 1-3 from range of TOF, we 
can construct the momentum of the recoil and its emis- 
sion angle with respect to the polarization direction. By 
fitting the angular distribution of the recoils for different 
momentum bins, we can extract the recoil asymmetry as 
a function of recoil momentum. 

The dependence on recoil momentum allows us to ex- 
tract information about the recoil terms while simulta- 
neously fitting for the tensor interaction. These different 
terms produce a different functional dependence of the 
recoil asymmetry on momentum. 

If we were to integrate over all momenta, the measured 
asymmetry Aspin as the spin polarization P is flipped 
can then be found from Eq. 1: 



_ w[e,p]-w[e,-p] 

^spin - w[9,P]+W[9,-P] ^ ' 

_ XiPArecoilC0s6> 

1 + cTx2 + cTx2Cos^6 

still ignoring recoil-order terms. An example of such a 
fit that could extract the quantity -P^rocoii is shown in 

Fig.m 

We wish to generalize this to a fit of the recoil asym- 
metry as a function of recoil momentum. In the absence 
of the Fermi function and of recoil-order corrections, we 
could define the recoil asymmetry as a function of recoil 
momentum ^spin [Pr] in terms of the kinematic functions 
of Appendix A: 
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FIG. 7: An example of the recoil asymmetry fit as a function 
of angle, with angle reconstructed from the impact location 
and the ion TOF. Here the asymmetry is for the charge state 
2 data, summed over all recoil momentum. 



^ W[0,P,Pr]^W[0,~P,Pr] 

^ w[e, p, Pr] + w[e, -P, Pr] 
UaAt - f7bT)Pcos9 



(9) 



fi - brfe - (a/3. + f )/2 



cr(/3 + /5COS2(0)) 



where all the fi functions depend on P^. However, we 
must properly include the Fermi function (which does not 
matter quantitatively) and the recoil order terms (which 
do matter). So instead we write the simple expression 
actually used for the fits as a function of angle, for each 
bin of recoil momentum: 




FIG. 8: (Color online) The dependence on recoil momentum 
of the experimental recoil asymmetries Ai. Charge states 1,2, 
and 3 are shown to be in statistical agreement. (The points 
for charge states 1 and 3 are offset horizontally for clarity.) 




after corn, 
no corr. 



2 3 4 5 6 

P ■, [MeV] 

' recoil 

FIG. 9: (Color online) Ai as a function of recoil momentum, 
determined with no correction for /3's and after the correction 
for /3^'s described in Section III CI The correction is small. 
Lines are drawn only to direct the eye. 



^spin [Pr'] 



PAi[Pr]cOS0 
l + cTF2[Pr]cOS^0 



(10) 



where F2 [Pr] now contains all the numerical integrations 
needed for the cos^(6') terms. These fits let us extract 
the experimental coefficient of cos(6') of the recoil asym- 
metry, >li[Pr]- We then will fit Ai[Pr] below, using the 
different dependence on recoil momentum of the new ten- 
sor physics terms and the recoil order terms, in order to 
extract the tensor physics terms At and bx- 

We show one experimental example of the spin-flip 
asymmetry ^spinl^*] as a function of cos{9) in Fig. [71 We 
make similar fits to ^spin as a function of binned recoil 
momentum from 0.5 to 5.0 MeV/c to extract Ai[Pr]. 



B. Correction for /3+'s striking MCP 

Since there is generally concern in precision /3 asymme- 
try experiments with the size of any corrections, we show 
in Figure |5| the results for Ai with and without the cor- 
rection from /3+ 's striking the electron MCP, as described 
in Section fll CI The final answer for the average ^recoil, 
shown below, becomes more negative by about 0.01 when 
the Z?"*" correction is made. We know this correction to 
better than 5% of its value. 



IV. RESULTS OF FITS TO THE 
EXPERIMENTAL Ai 



A. Experimental results for Ai[Pr] 

We show the results for Ai[Pr] broken down for the 
different charge states in Figure [8l The first three charge 
states can be seen to be statistically consistent. We there- 
fore simply take the weighted average over the charge 
states to consider the physics below. 



We now show results for fits to the extracted average 
v4i[Pr]. We try to take advantage of the different momen- 
tum dependence produced by the recoil-order corrections 
6m and d, and the non-SM tensor physics parameters At 
and bx- We consider different constraints on the param- 
eters, partly in hopes that in the future values of the 
recoil-order corrections may be available. 
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A. All parameters floating 



First we let all parameters float, 
points in the Ct + C'j, vs. Ct 



C'rj. plane, let &Af 
d float for each point, and compute the 



We consider various 
and 
We show 

the 90% CL limits in a contour plot in Fig. [TOl includ- 
ing constraints from other experiments. The best fit is 
for 5j\//(AA/GT)=-7.2±5.0, a negative value with rather 
large absolute value, and a small value of d/ A—Q±n . 
The square root of the reduced 'X^ /N is 1.37 for the best 
fit, so we expanded the error bars on the fit parameters 
appropriately. 



Tensor Limits (90%) 



u 
I 




-1.0 -0.5 0.0 0.5 

(c, + cy/c. 



FIG. 10: (Color) Exclusion plot at 90% confidence showing 
complementarity of the present constraints to other measure- 
ments. Allowed regions are inside: ''He P-v 0, |3|, red con- 
centric circles (hashed in between); polarization in ^*0, 
[3] (assuming no scalar interaction), blue solid vertical 
lines (hashed in between); Present work, letting 6m, rf, and 
both At and br float, black dashed hyperbolae. For br set 
to (consistent with J,]), weak magnetism bM / [AMgt)= 4.7 
± 4.7, and d left floating, then the green rectangular area 
indicates the present limit |(Ct — C^)/Ca| < 0.36. 



As a tool to discuss this result, we show the individual 
contributions from 6m, d, At, and 6^ in Fig- HH The 
momentum dependence of Ai from the non-SM tensor 
terms At and 6t is similar, so they are both made large 
with opposite signs to fit the data. The asymmetries 
produced by these non-SM tensor terms in this fit are 
therefore much larger in absolute magnitude than the 
experimental asymmetries. 

So a completely unconstrained fit does not produce 
competitive limits in a model with all chiralities of tensors 
possible and no constraints on the recoil order terms, 
because of the large number of similar degrees of freedom. 
It is also clear that the statistical precision of the data 
would be much better than the size of these error bars, 
if other constraints on the physics were applied. 




2 3 
P ., [MeV] 

recoil 

FIG. 11: (Color online) The dependence on recoil momen- 
tum of the recoil asymmetries for the best fl t with 6a/, d and 
both chiralities of tensor floating. yxV]V=1.37 for this fit. 
The cyan medium-dash line 'r.o.' is a small exact recoil-order 
correction given when recoil energy is included in the energy 
conservation equation. 



B. Constraining Ct + (7^=0 

We next use the relative positron polarimetry experi- 
ments comparing pure Fermi ^"^O and the pure Gamow- 
Teller branch 0] to imply that the Fierz interference 
term 6t cx: Ct + C't is very small. We continue to let 6m, 
d, and At float. 

The resulting constraints on Ct — C't can simply be 
seen by looking at the intersection of our limit hyper- 
bolae with the the Ct + C^=0 axis of Fig. [TOl A 
more detailed look at the minimum implies \{Ct — 



C'^)/Ca h0.42l^;^^ at 90% confidence, and the result- 
ing recoil order terms are bu / {AMct)=A%±.\\ and 
d/A—24:±4:l. Note here that 6m hes far outside the nu- 
cleon value. 




FIG. 12: (Color online) The dependence on recoil momentum 
of the recoil asymmetries for the best fit with Ct + C't=0 
(i.e. assuming the result of Ref. 0]), and 6a/, d and Ct — C't 
floating. The fit has a/xV^= 1-41, and the result is |(Ct — 
C't)ICa\ = 0.42155;;^^ at 90% confidence. 



The breakdown into the various contributions is shown 
in Fig[T21 Here it can be seen that 6m and At produce 
contributions to A\ with similar momentum dependence, 
so again it is difficult to fit them both simultaneously. It 
also is clear the the contribution to A\ from d is compar- 
atively constant with momentum, so d can be floated in 
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a meaningful fashion. 



C. Constraining Ct + C^=0 and constraining 6m 
from theory 

Here we keep Ct + C^=0 from the literature, and also 
constrain the weak magnetism term bM / {AMqt) from 
other physics. We assume it is given by the nucleon value 
4.7, with arbitrary error given by the full value 4.7. As we 
discussed in Section lTCl there are a number of cases in the 
literature where the other matrix element contributing to 
&M is not very large, and this range of 6jv/ easily covers 
all the cases in the literature. 

We list in Table H] the fit results for At and d for dif- 
ferent values of bu, in hopes that further knowledge of 
bM and d will eventually become available. We show the 
contributions to the fit from the different terms in Fig.lT^ 



0.015 ± 0.029 (statistical) ± 0.001 (systematic) ± 0.019 
(theory syst). The first systematic error is from the nu- 
clear polarization. The second systematic error is domi- 
nated by the error on bM / {AAdcr) of 100% of the nucleon 
value 4.7. 

Then, by using Eq. [51 we can interpret our result for 
Accoii to imply that CtC'j./C\ = 0.012 ± 0.022 ± 0.015. 
We show the resulting 90% confidence limits in Fig. [TO] as 
the green horizontal-hashed rectangular exclusion region 
showing \{Ct — C'j,)/Ca\ < 0.36 (with Ct + C'j, assumed 
0). It can be seen that under these assumptions, we can 
place constraints on tensor interactions coupling to non- 
standard Model right-handed neutrinos, i.e. on Ct — 
C'rp^ competitive with constraints derived from the fi-v 
correlation in ^He [1, 0]. 



6m/(AMgt) 


d/A 


At 




9.4 


-34 ± 42 


0.034 ± 0.032 


2.06 


4.7 


-24 ± 39 


0.015 ± 0.029 


1.90 





-14 ± 36 


-0.003 ± 0.027 


1.76 


-4.7 


-4 ± 34 


-0.022 ± 0.025 


1.64 


-9.4 


6 ± 32 


-0.040 ± 0.024 


1.56 


-14.1 


16 ± 31 


-0.059 ± 0.023 


1.51 


-18.8 


26 ± 31 


-0.078 ± 0.023 


1.50 


-23.5 


37 ± 31 


-0.096 ± 0.024 


1.52 


-18±11 


24 ± 41 


-0.074 ± 0.050 


1.60 



TABLE I: Results for fits to the dependence of the recoil 
asymmetry on recoil momentum for d/A and At, with N = 8. 
bM /AMqt is fixed at each of the values in the first column. 
We include the nucleon value of bM /{AMgt) = 4.7, and 
sweep through a large number of other values. The bottom 
line is the fit result if bM is allowed to float unconstrained; 
note that this value for bM is far outside the nucleon value. 



< 

E 
E 



on 



FIG. 13: (Color online) The dependence on recoil mo- 
mentum of the recoil asymmetries for the best fit with 
Ct + Ct~0 (i.e. assuming the result of Ref. 01), but 
bM/{AMGT)=4:.7±4.7, with d and Ct - C't floating. The 
result is At=0.015±0.029±0.019 (see text). 

We interpret the results of this table to mean At= 




D. Future improvements 



In order to improve the precision of this experiment, 
2nd-order recoil terms that might contribute to the 1+ 
to 2+ transition would have to be addressed. By adding 
efficient 7-ray detection to measure the recoils in coinci- 
dence with the 617 keV 7 ray, the experimental asym- 
metry for the excited-state transition could be measured 
separately. Then it could be included in the model for 
the total asymmetry. That would allow the extraction of 
new physics from the 1"^ to 0"^ transition free of the com- 
plications of the higher-order recoil effects of the excited 
state transition. This would be necessary to improve the 
sensitivity to below the 0.01 level. Other improvements 
would include implementing the better optical pumping 
used in Ref. [23j . which would improve accuracy by al- 
most a factor of 2 by making the polarization close to 
1.0. 

An accurate simultaneous (3 asymmetry experiment as 
a function of /3 momentum could also help constrain the 
recoil-order terms. 

Possible extensions of this experiment would include 
measuring the same quantity in ®^Rb [3l[ as a nuclear 
structure consistency test in the same shell. The ground 
state transition is an 82% branch, and the lower Q-value 
(3.4 vs. 4.7 MeV) also would help reduce dependence 
on recoil-order terms. The 3"*" to 2+ decay of the "^^K 
ground state could also be used to search for tensor 
interactions — although several 2nd-order recoil nuclear 
matrix elements are needed, calculations in the SD shell 
could be done with reliable estimates of the theoretical 
error. Recoil asymmetry measurements in ^Li would re- 
quire reconstruction of the momentum of the a particles 
emitted, and could contribute to searches for 2nd-class 
tensor interactions 32] along with fundamental tensor 
interactions. 
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V. CONCLUSIONS 

We have used atom trap technology to make the first 
measurements of the asymmetry of daughter nuclei with 
respect to the nuclear spin, as suggested by Treiman 
By measuring the momentum dependence of the asym- 
metry, we can constrain the recoil order induced tensor 
d independently. The similar momentum dependence of 
weak magnetism 5m and the non-SM tensor physics 6t 
and At produces correlations in their extraction. So in 
Section IIVI above we have considered the constraints on 
non-SM tensor interactions from the present experiment, 
using different assumptions from other experiments and 
from theory. 

If we make no assumption about tensor parameters 
from other experiments, and do not constrain bM and 
d, then the 90% CL constraints are shown in Fig. [TOl 
These are not competitive for tensors coupling to SM left- 
handed neutrinos, but their intersection with the vertical 
(Ct + C^) = axis shows the potential sensitivity of the 
technique. 

We believe our most useful result is to assume that 
tensors coupling to SM left-handed neutrinos do not ex- 
ist (i.e. Ct + Ct = 0) in accordance with Ref. [3], and 
let d float while constraining Bm /{AMot) to the nucleon 
value with 100% error at one sigma. Then we extract At 
= 0.015 ± 0.029 (stat) ± 0.019 (syst). Eq. Othcn im- 



plies CtC'j,/C\ = 0.012 ± 0.022 ± 0.015. The resulting 
constraints on tensor couplings to right-handed neutri- 
nos |(Ct — Ct)/Ca\ < 0.36 are also indicated in Fig. [TOl 
These constraints are complementary to those from the 
^He (i-iy correlation. We include Table I in hopes that 
future improvements in the knowledge of the recoil-order 
terms can be included. 

So our results place constraints on tensor interactions 
complementary to those from other experiments. The 
systematic error is dominated by uncertainty in extract- 
ing recoil order corrections from our data, and there are 
no serious experimental systematics. The statistical er- 
ror could be made considerably smaller with more count- 
ing and modest experimental improvements, which would 
cut the overall error by more than a factor of two even 
without theoretical guidance, and possibly by the order 
of magnitude that would make it complementary to the 
best measurements in nuclear, neutron, and pion decay. 
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APPENDIX A 

We present here the analytic results for the recoil spin 
asymmetry as a function of recoil momentum, assuming 
no Fermi function. These are obtained from the expres- 
sions for the full angular distribution in Refs. and [Toj 
by integration over v and (3 momenta. See Ref. (33j for 
more details. The units set mp—l. 

The angular distribution of daughter nuclei with re- 
spect to the nuclear spin as a function of their momentum 
Pr is 



W[Pr,0]dPrd{cOS Or) 



= [hiPr) + hrkiPr) - (a/3. + -j)h{Pr) + cT h{Pr) 

-PATf4{Pr)cOs{9r) + cT h{Pr)cOs'' {9r) 
+PbTf7{Pr)cOs{0r)]dPrd{cOS Or)- 

Integration of these kinematic functions over recoil mo- 



mentum Pr produces the terms in Eq. 1. We also in- 
clude here the normalization effect of the Fierz interfer- 
ence term 6^ and a similar effect Jdt in the asymmetry 
term, each multiplied by mp/Efj before integration over 
(3 energy. These terms are explicitly neglected in Eq. 1 
and in Ref. 

Note that the effect of the normalization term scaling 
/e is neglible in this work, because in Gamow- Teller de- 
cays Ai.ocoii=0 in the absence of tensor terms and the ef- 
fects of /e enter in higher order in the small tensor param- 
eters. That is no longer true for mixed Fermi/ Gamow- 
Teller transitions. 

The kinematic functions of momentum are given by 



h{Pr) 

( 

h{Pr) 

( 

h{Pr) 
h{Pr) 
MPr) 
h{Pr) 
h{Pr) 



{Pr - E^Pr + PfViSE^ + Pr' + P^ + E^{i - APr)) 



12{EI - Pr2)3 



{Pr - ElPr + P^fi^Et + P^hP^ - 1) - gg(3 + 8P,2)) 

12{EI - Pr2)3 
P^{l-El+P^f 
12{EI - Pr^Y 
EoP^{2 + El - P2)(i _ + P2)2 
6{Ei~P?r 
P^{2 + El - P2)(l ^El+ P^f 
Q{El~-P^f 
ElP^ - P^Eg + 2^oP2 
2{El-P^) 

{El + P2 + l){El P2 - 1) - ElPr{El - Pr') 

2{El-P^) 



where Eq is the maximum total energy of the /3, and 
where energies and momenta are in units of m/j. Inte- 
grating over the daughter nucleus momentum produces 
the terms needed in Eq. 1: 



Xi 



X2 = 



5{El ~ 6El + 3Eo + ^ + UEolnEa) 



4y^El^{2E^ - 9El - 8) + 
/B2^(4P4 _ 28EI - 81) + 15(6 + 



A, 



Ain = ln(£;o 



El 1) 



Note the typographical error in the equation for xi in 
Ref. 0. 



